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Whereas atom-molecule collisions have been studied with complete quantum 
^ ■ state resolution, interactions between two state-selected molecules have proven 

much harder to probe. Here, we report the measurement of state-resolved in- 

• i-H ■ 

^ ■ elastic scattering cross sections for collisions between two open-shell molecules 

that are both prepared in a single quantum state. Stark-decelerated OH radi- 
cals were scattered with hexapole-focused NO radicals in a crossed beam con- 
figuration. Rotationally and spin-orbit inelastic scattering cross sections were 
measured on an absolute scale for collision energies between 70 and 300 cm 1 . 
These cross sections show fair agreement with quantum coupled-channels cal- 
culations using a set of coupled model potential energy surfaces based on ab 
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initio calculations for the long-range non-adiabatic interactions and a simplis- 
tic short-range interaction. This comparison reveals the crucial role of electro- 
static forces in complex molecular collision processes. 

Rotationally inelastic scattering is one of the key processes underlying the exchange of 
energy between molecules (7, 2). In bulk systems, rotational energy transfer (RET) is respon- 
sible for the thermalization of state populations following a chemical reaction. In the dilute 
interstellar medium, inelastic collisions contribute to the formation of non-thermal population 
distributions that result in, for instance, interstellar masers (3). Accurate state-to- state inelastic 
scattering cross sections are essential ingredients for reliable models of chemical processes in 
combustion physics, atmospheric science, and astrochemistry. 

In molecular beam collision experiments, the ability to prepare molecules in a single rota- 
tional (sub)level prior to the collision using electric, magnetic, or optical fields has been im- 
perative to unravel the underlying mechanisms of molecular energy transfer. This has made 
scattering experiments possible at the full state-to- state level, and has resulted in the discovery 
of propensity rules for inelastic scattering (4), the stereodynamics of molecular collisions (5, <5), 
and quantum interference effects (7-9). The latest beam deceleration and acceleration meth- 
ods (10, 11) allow for the precise variation of the collision energy, resulting in the observation 
of quantum threshold effects in the state-to- state cross sections (72, 13). This wealth of stud- 
ies has contributed enormously to our present understanding of how intermolecular potentials 
govern molecular collision dynamics. 

Thus far these methods have mostly been used to study collisions of state- selected molecules 
with rare gas atoms. Yet, in most natural environments molecule-molecule interactions play a 
major role. For instance, space telescope observations of cometary water may reveal the pos- 
sible origin of water on Earth, but a conclusive interpretation requires accurate knowledge of 
RET in water- water collisions (14). Whereas atom-molecule scattering cross sections can now 
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be calculated routinely in excellent agreement with experiment (13, 75), much less is known 
about RET in molecule-molecule collisions (16). As opposed to an atomic target, a molecular 
scattering partner possesses internal degrees of freedom of its own, adding a level of complexity 
that can easily render ab initio quantum scattering calculations extremely challenging - if not 
impossible. This is particularly true for collisions involving radical species that are governed 
by multiple Born-Oppenheimer (BO) potential energy surfaces (PESs) with non-adiabatic cou- 
plings between them. Experimental data on bimolecular state-to- state cross sections is generally 
lacking, and kinetic models often use collision rate coefficients that are expected to be inaccu- 
rate (77). 

The study of molecule-molecule collisions at the ultimate quantum level has been a quest in 
molecular beam physics since it was established in the 1950s (78). Major obstacles exist that 
have prevented studies of state-to- state bimolecular scattering (79). The main challenge is the 
need for reagent beams with sufficient quantum state purity at the densities necessary to observe 
population transfer in one, or both, reagent beam(s). Thus far, experiments of this kind have 
only been possible using cryogenically cooled H 2 molecules as a target beam (20, 21). 

Here, we report the successful measurement of state-resolved inelastic scattering between 
two state- selected molecular beams. We have chosen the OH (X 2 II) + NO (X 2 II) system (22) 
for this purpose, as both open shell radical species are benchmark systems for the scattering 
of state- selected molecules with rare gas atoms that involve two BO PESs (23). Collisions 
between OH and NO involve eight interacting PESs, representing the full complexity of bi- 
molecular inelastic collisions (24). The OH-NO system serves also as a prototypical example 
of radical-radical reactions of fundamental importance in gas-phase chemical kinetics (25). We 
used a Stark-decelerator and a hexapole state selector in a crossed molecular beam configuration 
to produce reagent beams of OH and NO radicals with an almost perfect quantum state purity. 
The collision energy was varied between 70 and 300 cm -1 by tuning the velocity of the OH rad- 
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icals prior to the collision using the Stark decelerator, revealing the quantum threshold behavior 
of the state-to- state inelastic scattering cross sections. The unusually well defined distributions 
of reagent molecules allowed us to determine absolute scattering cross sections, which can nor- 
mally be determined only on a relative scale in crossed beam experiments. These cross sections 
showed fair agreement with a theoretical model for inelastic collisions between two 2 II radical 
species, based solely on an accurate description of the full rotational and open-shell structure 
of both radical species and their long-range non-adiabatic electrostatic interactions. This study 
revealed that inelastic scattering predominantly occurs at large intermolecular distances, even 
for the relatively high collision energies probed here. 

The crossed molecular beam apparatus used to study inelastic collisions between OH and 
NO radicals is schematically shown in Fig. 1 (26). A packet of OH radicals [X 2 U 3 / 2 ,v = 
0, j = 3/2, / (22), referred to hereafter as F 1 (3/2f)] with a tunable velocity in the 200 to 750 
m/s range was produced using a 2.6-meter long Stark decelerator (27). The velocity of the OH 
radicals was tuned by applying a burst of high voltage pulses to the electric field electrodes at 
the appropriate times. The state purity of the OH packets was such that less than 0.01 % of the 
OH radicals populated a lower A-doublet component of any rotational level. 

A beam of NO radicals with a fixed velocity was produced by seeding NO in a xenon 
carrier gas, and was passed through a 30 cm long electrostatic hexapole. NO radicals in the 
low-field-seeking Fi(l/2/) state were focused into the collision region, while molecules in 
the high-field-seeking Fi(l/2e) state were deflected from the beam axis. A 2 mm diameter 
beamstop and diaphragm were installed in the center of the hexapole and 10 mm downstream 
from the hexapole, respectively, effectively filtering out the Xe atoms from the molecular beam 
pulse. The resulting state purity of the transmitted NO F 1 (l/2f) beam was better than 99.0 %. 

The reagent beams of OH and NO were detected state- selectively in the collision region 
by a laser induced fluorescence (LIF) detection scheme. From calibrated LIF measurements, 
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the peak densities of the reagent packets were determined to be (2 ± 0.8) x 10 8 and (9 ± 3) x 
10 10 molecules/cm 3 for OH and NO, respectively. The collision induced populations in the 
Fi(3/2e), Fi(5/2e), Fi(7/2e), and the F 2 (l/2e) levels of the OH radical were measured at 
the time when both beams maximally overlapped in the beam crossing area. Depending on the 
inelastic channel, only a fraction of 10~ 4 to 10~ 6 of the OH radicals were inelastically scattered. 
Only final states of e symmetry were probed, as the Stark decelerator did not eliminate the initial 
population in the / states sufficiently. Detection of collision induced population in the Fi (3 / 2e) 
level was only possible by spectroscopically separating the magnetic dipole allowed transitions 
that originated from the Fi(3/2/) state (26). The insufficiently perfect state purity of the NO 
radical beam prevented the measurement of population transfer in NO. 

The collision signals were measured as a function of the collision energy, from which the 
excitation functions of the state-to- state inelastic scattering cross sections (shown in Fig. 2) were 
determined (28). The extremely well defined spatial distributions of the OH and NO packets 
allowed us to determine the scattering cross sections on an absolute scale. A cross section of 
90 ± 38 A 2 was determined for the Fi(3/2e) channel at a collision energy of 220 cm -1 , from 
which the absolute cross sections for all scattering channels and all collision energies were 
derived (26). The experimental uncertainty was limited only by the uncertainty in the measured 
value for the peak density of the NO packet. 

We found that collisions that populate the Fi(3/2e) level are most likely, and the cross sec- 
tion for this transition accounts for about 90 % of the total inelastic scattering cross section. 
The cross sections to populate the F 1 (5/2e), the F 1 (7/2e), and the F 2 (l/2e) levels show a 
clear threshold behavior; the collision energies at which these channels become energetically 
possible are indicated by vertical arrows. These cross sections show large qualitative differ- 
ences compared to the scattering of OH with atomic targets. The most striking difference is 
found in the relative contributions of the i<\(3/2e) and Fi(5/2e) channels to the total inelastic 
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scattering cross section. The role of the Fi(5/2e) channel, which dominates RET for OH-He 
and OH-Ne, gradually diminishes in favor of the F 1 (3/2e) channel in the series of targets He, 
Ne, Ar, Kr, and Xe (29). This behavior could be rationalized from the increasing well depth, 
anisotropy, and head-tail asymmetry of the two BO PESs (29, 30). The overwhelming domi- 
nance of the Fi(3/2/) —> F 1 (3/2e) quenching channel observed here for OH-NO reflects how 
a dipolar open- shell molecular scattering partner rather than a spherical atomic partner governs 
the collision dynamics. 

In order to interpret our experimental results, we constructed a model for the scattering of 
two molecules in an open shell 2 II state. In contrast to scattering of OH or NO with rare gas 
atoms, ab initio calculations of multiple anisotropic PESs with their non-adiabatic couplings for 
OH-NO are beyond the capabilities of current theoretical methods. Coupling of the S = 1/2 
electron spins gives rise to singlet (S = 0) and triplet (S = 1) potentials, which describe differ- 
ent short-range exchange interactions. There are four spatially distinct electronic states for each 
spin state, which are degenerate at long range and for linear geometries and which are coupled 
by non-adiabatic interactions. Nuclear derivative couplings with respect to all nuclear degrees 
of freedom exist between these states. Ab initio studies of the OH-NO complex (31) focused 
on the region where the chemical reaction OH + NO —> H + N0 2 takes place, but considered 
only the lowest adiabatic potential for the singlet state. Even if we were able to compute all the 
relevant adiabatic PESs, there would be no simple recipe to take the non-adiabatic couplings 
between the PESs into account. 

In our model we exploit the hypothesis that the processes with the largest cross sections are 
governed by couplings that occur at relatively large OH-NO separations, beyond the HONO 
well region. As opposed to the short range interactions, the long-range parts of the PESs can be 
calculated accurately by ab initio methods. We neglected the complicated short-range behavior 
of the PESs and replaced it with an isotropic repulsion term. However, we accurately calcu- 
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lated the long-range PESs that are governed by first-order electrostatic interactions between the 
dipole, quadrupole, and octupole moments of the collision partners (26). Moreover, we included 
isotropic dispersion and induction terms. The intermolecular Hamiltonian contained the usual 
radial and centrifugal kinetic energy operators, and the full 4x4 matrix of diabatic interaction 
potentials (26, 32). Due to the non-cylindrical symmetry of the 2 II ground states of both the 
OH and NO radical, the off-diagonal elements of this matrix provided by the quadrupole and 
octupole moments of both radicals contain important couplings between the 2 n 3 / 2 and 2 Ui/ 2 
states of both species. 

The cross sections that were obtained from these model PESs by full coupled-channels 
calculations are shown as solid lines in Figure 2. Fair agreement between experiment and 
theory was obtained, in particular considering the simplistic approximations for the short-range 
PESs that were made. The absolute value for the cross section of the dominating Fi(3/2e) 
channel, as well as the relative strengths of the inelastic channels, are reproduced well by the 
model calculations. The cross section for the spin-orbit changing F 2 (l/2e) channel, as well as 
the cross section for the F 1 (3/2e) channel at low collision energies, are overestimated by the 
model. 

We tested the sensitivity of the model calculations with respect to changes in the short-range 
repulsion term (26). We observed that the cross section for the parity-changing F 1 (3/2e) chan- 
nel is governed exclusively by the long-range electrostatic interaction; its value is converged 
to within a few percent. About half of this large quenching cross section originates from colli- 
sions with impact parameters exceeding 12 a . The F 1 (5/2e) channel is also mainly determined 
by the long-range forces, although its cross section varies by 10 to 25 % upon changes in the 
short-range model parameters (26). The weak F 1 (7/2e) and F 2 (l/2e) channels show larger 
variations, and more realistic short-range PESs are required to accurately predict their cross 
sections. 
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Our model also predicts the final states of the NO radical that are populated in coincidence 
with RET in the OH radical, but that cannot be probed with the present experimental arrange- 
ment. The dominant Fi(3/2e) quenching transition in OH is accompanied by the inelastic 
channels in NO as given in Table 1 for various collision energies. The general scattering be- 
havior can be understood from the terms that lead to inversion parity changing or conserving 
collisions with respect to both collision partners. The OH-NO dipole-dipole interaction re- 
sults in transitions that either change or conserve inversion parity in both OH and NO, whereas 
the dipole-quadrupole and quadrupole-dipole terms also allow for an inversion parity changing 
transition in OH or NO only. The largest cross sections that accompany the i<\(3/2e) channel 
in OH are found for the dipole-dipole dominated F 1 (3/2f) and F 1 (l/2e) channels of NO, and 
for the dipole-quadrupole dominated Fi(l/2/), Fi(3/2e), and F 1 (5/2f) NO channels (33). 

Our experiments show that the main mechanisms of RET in the reactive, chemically com- 
plex system studied here are captured using a model for the long-range interactions alone, pro- 
vided that the full monomer Hamiltonians and all relevant long-range non-adiabatic couplings 
are taken into account. Even at relatively high collision energies, the inelastic scattering events 
with the largest cross sections predominantly occur at large intermolecular distances where the 
interaction potentials can be calculated accurately. The success attained here implies that re- 
liable predictions for state-to-state scattering cross sections can now be made more generally 
for complex molecular systems involving radicals, helping to solve urgent scientific questions 
in, for instance, astrochemistry. Ultimately, new electronic structure methods that include the 
chemically reactive short-range potentials and non-adiabatic couplings are required to elucidate 
the exact mechanisms of radical-radical collisions. 
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Table 1: Predicted state-to- state cross sections (A 2 ) for RET in NO (initial state F 1 (l/2f)) that 
occurs in coincidence with the F 1 (3/2f) —> F 1 (3/2e) transition in the OH radical at collision 
energies of 10, 20, 50, 100, 200, and 300 cm -1 . The inversion parity of the molecular levels 
is given. For comparison, the elastic cross section corresponding to the OH (F 1 (3/2f) —> 
Fi(3/2/))-NO(Fi(l/2/) -> Fi(l/2/)) channel is given in the last row of the table. 
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Figure 1 : Schematic representation of the experimental setup and the energy level schemes of 
the OH and NO radicals. A state- selected and velocity tunable beam of OH radicals produced 
using a 2.6-meter-long Stark decelerator was crossed with a hexapole state- selected beam of NO 
radicals. Both radical species were detected state- selectively using laser induced fluorescence 
(LIF), with total fluorescence intensity measured using a photomultiplier tube (PMT), and the 
spatial distribution of both reagent molecular packets imaged onto a charge-coupled device 
camera. Typical images of the OH and NO packets are shown in the upper and lower insets, 
respectively. The mean speed of the OH radical packet was precisely known from the settings 
of the Stark decelerator. The collision energy was calibrated from the NO beam speed measured 
via a second LIF detection zone located 30 cm downstream from the collision area. The X 2 II 
electronic ground states of the OH and NO radicals are split into two rotational manifolds due to 
the spin-orbit interaction. The manifolds with lowest energy [\Q\ = 3/2 for OH and \Q\ = 1/2 
for NO] are labeled F\. The energy splittings shown between the A-doublet components of each 
rotational level are greatly exaggerated for clarity. 
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Figure 2: Comparison of the collision energy dependence of the measured (data points with 
error bars) and calculated (solid curves) state-to- state inelastic scattering cross sections of OH 
Fi(3/2/) radicals in collision with NO i<\(l/2/) radicals. The cross sections were first mea- 
sured relatively with respect to each other, and the vertical error bars indicate combined esti- 
mates of both statistical and systematic errors (2a) (26). The vertical axis was then put on an 
absolute scale by a measurement of the absolute cross section for the Fi(3/2e) channel at a 
collision energy of 220 cm -1 . The cross sections were computed on an energy grid of 10, 20, 
30, 40,..., 320 cm -1 . The cross section for the dominant Fi(3/2e) channel converged to within 
a few percent; the cross sections for the weaker channels vary by 20 to 50 % depending on 
changes in the short-range part of the theoretical model (26). 
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